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Abstract 

t-H ■ 

" At high values of the pair invariant mass the differential cross section for top-quark 

pair production at hadron colliders factorizes into soft, hard, and fragmentation func- 
tions. In this paper we calculate the next-to-next-to-leading-order (NNLO) corrections 
to the soft function appearing in this factorization formula, thus providing the final piece 
j — ■ needed to evaluate at NNLO the differential cross section in the virtual plus soft approx- 

imation in the large invariant-mass limit. Lechnically, this amounts to evaluating the 
vacuum expectation value of a soft Wilson loop operator built out of light-like Wilson 
lines for each of the four partons participating in the hard scattering process, with a 
certain constraint on the total energy of the soft radiation. Our result turns out to be 
surprisingly simple, because in the sum of all graphs the three and four parton con- 
tributions multiply color structures whose coefficients are governed by the non-abelian 
exponentiation theorem. 



1 Introduction 



The pair invariant mass distribution is an important observable in top-quark pair production 
at hadron colliders. Especially interesting is the high invariant-mass region of the distribution, 
which could potentially be distorted by new physics without disturbing the good agreement 
between Standard Model and experiment for the total cross section. The phenomenological 
importance of this distribution motivates special attention to its calculation in the Standard 
Model. 

Theory predictions of the invariant mass distribution typically rely on the next-to-leading 
order (NLO) computations of two-particle inclusive cross sections carried out in [1J, supple- 
mented with soft gluon resummation at the level of next-to-leading-logarithmic (NLL) [2JE] or 
more recently next-to-next-to-leading-logarithmic (NNLL) [I] accuracy. The resummed cal- 
culations take into account higher-order logarithmic plus-distribution corrections related to 
gluon emission in the soft limit z = M 2 /s — > 1, with M the pair invariant mass and \/~§ the 
partonic center-of-mass energy. However, except at restrictively high values of the invariant 
mass where r = M 2 /s — > 1, with yfs the hadronic center-of-mass energy, the dominance of 
these logarithmic corrections over 5(1 — z) corrections relies on the mechanism of dynamical 
threshold enhancement studied in [5HS]. In fact, numerical comparisons of exact NLO results 
with the leading terms in the soft limit show good agreement only if the delta-function terms 
are included to achieve a full virtual plus soft approximation [9j. Starting at NNLO, such a 
virtual plus soft approximation is not achieved through expansions of NNLL resummation for- 
mulas alone, and even though the delta-function contributions are formally N 3 LL corrections 
the only way to know their size for certain is to calculate themQ Moreover, the soft gluon 
resummation mentioned above uses the generic counting m t ~ M for the top-quark mass, 
whereas at truly high values of the invariant mass the counting m t <M should be used, and 
resummation must also take into account logarithms of the ratio m t /M. 

In a recent paper [12], we have set up a factorization formalism appropriate for the 
invariant-mass distribution in the simultaneous soft and small- mass (m t <C M) limit. Schemat- 
ically, the factorization is of the form 

da = Tr [H S] ® D <g> D + 0(1 - z) + O Qj) . (1) 

The hard function H and the soft function S are matrices in the space of color-singlet operators 
for (qq, gg) — > ti scattering, evaluated with m t = 0, while the D are perturbative heavy-quark 
fragmentation functions containing the dependence on m t . Given this form of factorization, 
one can derive and solve renormalization-group (RG) equations for the individual functions to 
resum soft logarithms as well as those depending on m t /M. Just as importantly, the individual 
functions appearing in (Tj[|) are much easier to calculate in fixed-order perturbation theory than 
the hard and soft functions needed for soft gluon resummation for generic m t . In fact, the 
fragmentation function is completely known to NNLO accuracy [13] . and the higher-order 
virtual corrections to two-to-two scattering needed to extract the contribution of the NNLO 

1 The delta- function terms and also non-singular terms in the soft limit would be included in a full NNLO 
calculation, but current NNLO results are limited to certain qq-initiated contributions to the total cross 
section [TUlfTT] . 
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hard function to the differential cross section are also known [TIHTS] . The only missing piece 
needed to obtain at NNLO a full soft plus virtual approximation in the limit of large invariant 
mass is the soft function S. The calculation of this soft function to NNLO is the subject of 
this paper. 

Our main motivation for this calculation is its eventual impact on the phenomenology of 
top-quark pair production. However, we consider it an interesting problem even apart from 
this. The soft function we deal with here is related to double real emission corrections to 
massless two-to-two scattering, and at the technical level is defined as the vacuum expectation 
value of a Wilson loop built out of four light-like Wilson lines. While a number of soft functions 
have been calculated at NNLO in the literature, these all involve either two [T§l - f2"5] or three [2E] 
Wilson lines. Compared to those cases, the four- Wilson-line soft function is characterized by 
the added conceptual complication of a non-trivial matrix structure in color space. One might 
also expect it to be computationally much more complicated because, unlike the case studied 
in [26], graphs with attachments of gluons to three Wilson lines do not vanish and in general 
are complicated functions of two non-trivial scalar products. However, we find that in the sum 
of all diagrams such three-parton contributions multiply a color structure whose coefficient 
is determined by the non-abelian exponentiation theorem [2T], l2"8]. In particular, the bare 
function does not contain a three parton term with the antisymmetric color structure of a 
three-gluon vertex. This is an expected result for the IR divergent pieces of the bare function, 
as a consequence of the form of the RG equation derived in [12]. In particular, the momentum 
dependence in the anomalous dimension governing this RG equation is inherited from the 
anomalous dimension for scattering amplitudes and is of the dipole type at least to NNLO, 
which follows from results in [29], [30] (and may even be true to all orders, as conjectured 
in [3TH33] ). For the IR finite pieces it is perhaps slightly unexpected that the three parton 
terms are determined by non-abelian exponentiation, but this is nonetheless the case due to 
cancellations among certain diagrams. 

The remainder of this paper is organized as follows. First, in Section [2] we give the precise 
definition of the soft function calculated in this paper and we also review the NLO calculation 
from [12] , using it as a means of introducing some of the formalism related to the color-space 
matrix structure. Then, in Section[3]we present an expression for the bare NNLO soft function 
as a sum over legs and give explicit results for the component integrals and matrix structures 
appearing in this sum. We also describe cross checks, both with the two- Wilson-line integrals 
calculated in [23] and with the non-abelian exponentiation theorem. Finally, in Section HJ we 
discuss the renormalization procedure and explain how this provides a further cross check on 
our result. We conclude in Section El relegating several details of the calculation along with 
the final results for the NNLO soft function to the appendix. 

2 Definitions and the NLO calculation 

We define the soft function for the pair-invariant mass distribution as in [3], adapted to the 
massless case by replacing time-like velocity vectors by light-like ones [12]. The basic objects 
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for the soft function are firstly the Wilson lines 



Si(x) = V exp (ig s j dsrii- A a {x + sm) T^j , (2) 

where the V refers to path ordering, and secondly the Wilson loop built out of these objects, 

O s {x) = [iS' ni /S l n2 /S l ra3 iS' n4 ] (x) . (3) 

We have used the notation of [3U|35], where the bold-face indicates that the objects T t a are 
matrices acting on color structures specific to the type of partons participating in the two-to- 
two scattering process. This notation allows us to describe simultaneously the case where the 
Wilson lines Si associated with the four partons are in any combination of the fundamental (for 
quarks) or adjoint (for gluons) representations. In this paper we have in mind applications to 
top-quark pair production in the soft limit, and so will give results appropriate for (qq, gg) — > ti 
scattering. 

The Wilson-loop operator ([3]) takes into account the coupling of soft gluons to the external 
partons within the eikonal approximation. The soft function is related to the contribution of 
these gluon emissions at the level of the squared amplitude. In [12J, it was defined through 
the Fourier transform of a position-space soft function evaluated in the parton center-of-mass 
frame. In the present work, we will find it more convenient to work directly with the following 
momentum-space representation! 

S{uM/M 2 ^) = -L^ ( |O](0)|X s > (X s \O s (0)\0) 5(u - {n x + n 2 )-p Xa ) , (4) 
dR x a 

where cIr = N in the qq channel, and d,R = N 2 — 1 in the gg channel, with N the number of 
colors. As usual, X s refers to a final state built up of any number of unobserved soft gluons. 
It is clear from (J4J) that the soft function depends on the single dimensionful parameter u 
(and /i, upon renormalization) , as well as the scalar products rii ■ nj. Although the results 
we give later on can be used to construct the soft function for arbitrary velocity vectors, we 
have defined it in the natural way for two-to-two scattering. In that case there are only two 
independent scalar products and thus one non-trivial dimensionless ratio, which we have chosen 
as rti -n 3 /ni-n 2 = n 2 - n^/rti -n 2 = —ti/M 2 . Our notation is then in direct correspondence with 
the Mandelstam variables used in [12]. It further implies that rti ■ n 4 /?T, 1 • n 2 — n 2 ■ n 3 /rii ■ n 2 = 
1 + ti/M 2 = —U\/M 2 and ni ■ n 2 = n 3 ■ n 4 = 2. In the parton center-of-mass-frame, the delta 
function constrains the energy of the soft radiation to 2E S = u, and it is an easy matter to 
show the correspondence with the position-space definition used in [12J. 

In order to study the higher-order corrections which are the subject of this paper, we define 
expansion coefficients of the bare soft function in d = 4 — 2e dimensions as 

s b . re =s<°> + sg. + 2 sSL + ■ ■ • • (5) 




2 Note that this definition differs from [T^] by a factor of y/I, which we have chosen to omit here. The 
Laplace-transformed function in (jJTJ), on the other hand, coincides with the one introduced in |12j . 
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In the above equation we have expressed the bare coupling constant in terms of the 
renormalized one in the MS scheme: the relation between the two is Z a a s n 2e = e~ e7B (47r) e a;s 
with Z a = l—f3 a s / (47re) and (3 = 11/3N— 2/3rty, with rii the number of light flavors. The soft 
function for massless partons was obtained to NLO in [12]. We end this section by reviewing 
the elements that go into that calculation. This gives us an opportunity to introduce the 
aspects of the color-space formalism needed in this work, and to extend our previous results 
to the depth of the e-expansion needed for the NNLO calculation. 

The color-space formalism provides a means of representing soft gluon emissions from 
external quarks and gluons in a unified way. The important point for the case at hand is that 
these soft emissions can mix the possible color-singlet structures appearing in the two-to-two 
scattering amplitudes. To explain the matrix structure relevant for this mixing, we must first 
define a color basis for (q ai q a2 , g ai g a2 ) — > t a H ai scattering, where the {a} label the color indices 
of the partons with velocity n^. We work in the s-channel singlet-octet basis 

V°l ){ a } ~ u a 1 a 2 ^a 3 a A , \<-2 ) { a } ~ i a 2 a 1 i a 3 a 4 J 

( r 99\ _ XO.ia 2 s. (99\ _ • ea x a 2 c ,c f ag\ _ m x a 2 c ,c 

\ C 1 ){a} - ° a 3 a i ' \ C 2 ) {a } - l J r a 3 a 4 > l C 3 ) { a } ~ « l a 3 a 4 ■ \P) 

We view these structures as basis vectors \ci) in the space of color-singlet amplitudes. Inner 
products in this space are defined through a summation over color indices as 

( c | c > = y( c \* ( c ) . ( 7 ) 

{a} 

This inner product is proportional but not equal to 5u, so the basis vectors are orthogonal 
but not orthonormal. The soft function matrix elements are defined as 

S IJ =(c I \S\cj) . (8) 

The soft function is thus a two-by-two matrix in the qq channel, and a three-by-three matrix 
in the gg channel. At leading order (LO), a short calculation yields the result 
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(9) 



At NLO and beyond gluon emissions from eikonal lines are associated with factors of TJ. One 
determines the matrix elements (IS]) by using a set of rules which dictate how the act on the 
color indices {a} of the color structures (J2J) and performing the sum over colors indicated in 
d?]). These rules are as follows [34"]l35]: if the z-th parton is a final-state quark or an initial-state 
anti-quark we set (Z^ c )& a = £& a) for a final-state anti-quark or an initial-state quark we have 
{T^) ba = —t c ab) and for a gluon we use {T^) ha = if abc . In (@| there is also a conjugated operator 
O], which contains conjugated color matrices acting to the left. However, we can drop 
the daggers and let them act to the right, using the relation 

(M\T^\M'} = (I™ \M))^ \M') = \M) ] {T? \M')) = (M\T?\M') , (10) 
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Figure 1: Diagram contributing to the NLO soft function. 



with \M) and being two arbitrary vectors in color space. This can be easily understood 
by noting that T" 1 is a Hermitian operator. To show explicitly that this relation holds, it is 
sufficient to work out the case where i — 1 and \M) and are two basis vectors |oi, . . .) 
and . . .). We then have 



(7? h,...}^ |& 1; ...> 
l«i,.--> + (^i c |6i,...» 



((T l % iai \d 1 ,..»^b 1 ,..) = (T i y diai 8 bldl 
\a>i, ■ ■ •} ((7\ c )di&i \di, •••)) = (7T) dlfel 5 aidl 



Note also that if there is more than one T matrix acting on the same Wilson line, the conju- 
gation reverses their order, e.g., 



(2? 2? \M)) ] \M') = \M) ] (2*2? \M')) 



(11) 



This is crucial for the cancellations of certain contributions which we will encounter later on. 

The formalism is now in place to discuss the NLO calculation. This requires us to evaluate 
the type of diagram shown in Figure[T]and perform a sum over attachments to the different legs. 
The diagram represents the connection of gluons to Wilson lines for two distinct particles % 
and j. Diagrams involving virtual corrections or two attachments to the same line are scaleless 
and vanish in dimensional regularization. To evaluate the diagram we associate each gluon 
emission with the Feynman rules for eikonal attachments following from (T5]) and set the cut 
propagator (for which we use Feynman gauge) on shell with positive energy. This results in 
an expression proportional to the integral 



Ix(u},a,ij) = 
where [dk] = d d k 5{k 2 ) 6{k°) and 



[dk] ni • 5{U ~ U ? ■ k) = Tr^e— a," 1 -* 



(12) 



1 



TIq 71/ <i ' Tli j 

2n ■ rii n ■ rij 



(13) 



For top-quark production, where n| = {n\ + n 2 ) 2 = 4, we have ayi = 034 = 0, 013 = a 2 4 = 
—Ui/M 2 , a 14 = a 2 3 = —ti/M 2 . Later on we will describe cross-checks of our result with 
the simpler NNLO soft functions for Drell-Yan [19] and electroweak boson production at 



5 



large pt [26J. These involve the same two- Wilson-line integrals but with a = and a = 1 
respectively. 

The integral in ffT2]) can be parameterized in terms of the gluon energy and an angular 
integral which is of the type considered long ago in [36]. We have rederived the result using 
the light-cone coordinate decomposition from [22], which turns out to be especially convenient 
for the NNLO integrals considered later on. Either way, the result for the stripped integral is 

2e elE Y(— e) 

h{a) = T(l- 2e) (1 " ^ 2Fl( " e ' 1 " e ' a) ' (M) 
which can be expanded in e using 

2 F 1 (-e, -e, 1 - e, a) = 1 + # 2 (a)e 2 + (H 3 (a) - # 2 ,i(a)) e 3 

+ (H 4 (a) - H 3A (a) + if 2 ,i,i(a)) e 4 + • • ■ . (15) 

The functions if indicate Harmonic Polylogarithma^l (HPLs) [37J. Here and elsewhere in the 
paper we have used the Mathematica package HypExp [3B] in expanding the hypergeometric 
functions and manipulated the resulting HPLs with the package HPL |39j . 
We obtain the bare NLO soft function through the following sum over legs: 



legs 



t (u) 2e 7i(ai2) + w ™ 7i(ai2) + 2w ^ 7i(ais) + 2w ^ /i(au) ) • (16) 



We have taken into account the relations between the and also the explicit form of the color 
matrices given below to simplify the sum. The matrix structure in color space is obtained by 
evaluating the matrix elements of 

= — T 4 .T„ (17) 

as in OH]). Results were given in [I], and we reprint them here for convenience. In the qq 
channel they read 

(i) (i) C F (4N 2 
12 34 4iV -1 



CV (0 I 

2 



(1) (1) 

"-"•"111 2C,-« 



3 We write the HPLs in the compact notation which eliminates the zeros in the weight vector by adding at 
the same time one to the absolute value of the previous index to the right; for example H%(a) = H(0, l;o), 
Hi{a) = #(0,0, 0,1; a), #2,1,1(0) = #(0, 1, 1, 1; a), etc. 
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w {1) 



W 



(1) 
23 



Cp_ I —N' 
2N\-N 1 



while for the gg channel they are 



w 



(i) 

12 



(AN 2 \ 

ON 2 
\ N 2 — AJ 



w 



(i) 

31 











\ 









TV 2 -4 
47V 2 , 



W 



(1) 
13 



(1) 
21 



1 



/ o 



o \ 



AN 

AN N 2 N 2 -A 
\0 N 2 - A N 2 - A) 



w 



(i) 

14 



W 



(1) 
23 



/ 

-AN 

\ 



-4JV 
N 2 



-(N 2 - 4) 



\ 

-(N 2 -A) 
N 2 — A J 



18) 



(19) 



As usual, C F = (iV 2 - 1)/2N and = N. 

At NLO the renormalized function in the MS scheme can be obtained from the bare 
function simply by dropping the poles. More formally, we need to multiply it on both sides by 
a UV renormalization matrix. We describe the formal procedure in more detail in Section HJ 
after performing the NNLO calculation in the next section. 



3 The bare soft function at NNLO 

In this section we calculate the bare soft function at NNLO. We find that in the sum of all 
diagrams the result can be written in the form 

S bL = -{^) Yl ( Yl w if + w i% a lk ) + wf kl Igfaj, a kl ) 

legs \ra=2 



(20) 



2 (w^ + wity (/ 2 (a 12 ) + C A I 6 (a 12 ) + C A I 7A {a 12 ) + C A I 7>2 (a 12 )) 



= - (^Y e 

CO \U)J 

+ Aw[^ (l 2 (a 13 ) + C A I 6 (a 13 ) + C A I 7:1 (a 13 ) + C A I 7 , 2 (a 13 )) 
+ Aw { S (T 2 (a u ) + C A I 6 (a u ) + C A I 7>1 (a u ) + C A I 7>2 (a u )) 



+ 2 (w$ + u^?) I 3 {a 12 ) + Aw^ (/ 3 (oi 3 ) + h{a 13 )) + Awf} (J 3 (a 14 ) + / 4 (a 14 )) 




Figure 2: Two- Wilson-line integrals required in the calculation of the NNLO soft matrix. 
+ Aw { S h{a 12 ) + (w<$ + u4?) (14(012) + 2/5(012)) 

+ 2ti>g> (/ 4 (a 13 ) + 2/5(013)) + 2w$ (7 4 (a 14 ) + 2/ 5 (o 14 )) (21) 
+ 4 (w^ 3 + ioJ^J / 8 (ai2, 013) + 4 (logt + W324) /s(ai2, «i4) + 8w { Si h(a>i3, au) • 

Three types of basic diagrams contribute to this sum, depending on whether the gluons attach 
to two, three, or four distinct Wilson lines. We organize this section by discussing each type of 
diagram in term, and give explicit results for the color factors and integrals appearing in (|2T|) . 
The results for the three and four parton diagrams turn out to be surprisingly simple, because 

(8) 

the non-abelian exponentiation theorem constrains the coefficients of the color structures 
and wfh in fTSOj) . We discuss this further in Appendix [B] 

3.1 Two- Wilson-line integrals 

The subset of non-vanishing two- Wilson-line integrals is familiar from other calculations of soft 
functions to NNLO [THJEHEB] ■ The relevant Feynman diagrams are shown in Figure [2j Con- 
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verting these diagrams into integral expressions is straightforward, and leads to the following 
set: 



I 2 {u,aij] 



7T (C A (5-3e)-2^(l-e)) 



£(2-6) 
r(4 - 2e) 



rij ■ rij 5(u — no ■ k) 
rii ■ k rtj ■ k (k 2 ) 1+e 



[dk] [dl] 
[dk] [dl] 
[dk] [dl] 
[dk] [dl] 
Ii 7 i(u),a,ij) = — 7r 1_e 9R[e" 



(rij • rij) 2 5{oj — n ■ (k + I)) 
Hi ■ k rij ■ (k + I) rij ■ I rij ■ (k + I) ' 

{rii ■ rij) 2 8{lo — n ■ (k + I)) 

(rii ■ rij) 2 5(oj — n ■ (k + I)) 
rii ■ k rii ■ {k + I) rij ■ k rij ■ (k + 1) ' 

rii ■ n j n i ■ (I — k) 5{uj — n ■ {k + I)) 
rii-km-ik + l) rij ■ (k + I) (k + I) 2 

r 2 (i + e )r 3 (- e 



r(-2e) 



[dq] 



rii ' n j 



2rii ■ q rij ■ q 



l+e 



5(u -n -q) 



^7,2 (w,Oi_ 



[dk] [dl] 



rij rij ■ (k + 21) 8(co — n Q ■ {k + /)) 



rii ■ k rij ■ (k + I) rij ■ I (k + /) : 



(22) 



The prefactors in J 2 and I? t i arise from the internal loop integrals, see for instance the discus- 
sion in [19],[26]. To evaluate the phase-space integrals we use light-cone coordinate techniques 
and parameterizations described in the Appendix of [26]. It is then relatively straightforward 
to derive results for the integrals in terms of hypergeometric functions, with the exception 
of -Z5 and I-j^i which we evaluate as an expansion in e. In Appendix |A] we derive the results 
for the integral ^5 as an example of the calculational procedure. Defining stripped integrals 
according to (with n > 1) 



I n {w, dij / 



71 



2-2<e a -2e~f E -l-4e 



IjJ 



Inidij) 



the explicit results can be summarized as 
8e 2e lE r 2 (- e ) r(-2e) 



h(d) 



/ 5 (a) = (l-a) 



r(l - 2e) T(l - 4e) 

1 1 /7tt 2 



l-a)- 2e [ 2 Fi(-e,-e,l-e,a)] , 



-2c 



6 



62 

- 4H 2 (a) ) + — C 3 - 12H 3 (a) + 4ff 2)1 (a) 



(23) 

(24) 
(25) 



V 14tt 2 

i ' I — + — H 2 (a) - 36F 4 (a) - 4# 2 , 2 (a) + 12# 3 ,i(a) - 4# 2>lil (a; 



J 4 (a) 



(26) 
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2e 7 B rif—A r(-2A 

Ha) = r(2 _ 2e)r(1 _ 4e) (1 - a)- 2 ^ 2 F,(-2e, -2e, 1 - 2e, a) , (27) 
I 2 (a) = - (C A (5 - 36) - 2^(1 - 6)) ^(^t^r^^ Ha) , (28) 



7r,<«) = - ^^r(-^(-^(i + o (1 _ a) -, 2Fi( _ 2e> _ 2e> , _ e> a) 



(29) 



I 7 , 2 (a) = (1-a) 



-2c 



2 1 /^7r 2 \ 139 

+ - f — - 6tf 2 (a) J + — C 3 - 8tf 3 (a) + 12F 2jl (a) (30) 



V 23tt 2 

+ — H 2 (a) - 12# 4 (a) - 4# 2 , 2 (a) + 16# 3 ,i(a) - 24# 2)1 , 1 (a) ) 



To expand these in e we use (1151) along with 

2J Fi(-2e, -2e, 1 - e, a) = 1 + 4# 2 (a)e 2 + (AH 3 (a) - 12# 2>1 (a)) e 3 

+ (AH 4 (a) + 4# 2 , 2 (a) - 12H 3jl (a) + 36H 2jl>1 (a)) e 4 + • • • . (31) 

These results can be related to those for a position-space soft function calculated in [24J, and 
we have found full agreement with that work. They furthermore agree with results from [19] 
for a = and from [26J for a = 1 as two special cases. 

We must also evaluate the color factors. For these we find 



(2) (1) 

W ij = w u 



(3) _ -1- rpa rpb rpa rpb 

ij d R ' ' 3 1 ' 



v v dji 1 1 3 3 13 2 %3 

II 

1 1 ""'i ■ '/ ^ "■' " I ' '"J ■■) " './ 



-1? = «™ = ^ = i- *F bc T?T*T<=^ wff . (32) 



The relation between and along with the result that I 3 = 7 4 /2 — 7 5 , ensures that 
the bare function satisfies the non-abelian exponentiation theorem. We discuss this further 
in Appendix [Bl Results for the NLO matrices were given in (fIBl) and (fl9i) . In the qq 
annihilation channel the remaining matrices evaluate to 

(3) (3) C F 

w{ 2 ' = w& = — 



(3) (3) C F 



2 
10 



N 2 - 


'1 





AN 2 J 


N 2 -2 


N 2 ~2 
2N 
N 4 -3N 2 +3 


2N 


4AT2 



w 



(3) 
14 



W 



(3) 
23 




(33) 



whereas in the gluon fusion channel they are 
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(3) 
12 



10 



(3) 
31 
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(3) 
13 



2 
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/ 2iV 2 \ 
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TV 2 



V 

/ (iv 2 -i) 2 n 

' 2TV 2 U 

U 4TV 2 



(3) 
24 



Ca 
2 



/ 1 



TV 
4 

TV 2 -4 
. 47V 



o \ 



TV 



47V 4 . 

TV 
4 

7V 2 +2 
8 

7V 2 -4 



TV 2 



47V 
TV 2 — 4 



(TV 2 — 2) (TV 2 — 4) 
87V 2 



(3) 
™14 



(3) 
23 



Ca 
2 



_TV 
4 

TV 2 -4 
47V 



_7V 
4 
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3.2 Three- Wilson-line integrals 

The three- Wilson-line integrals are of two basic types: the abelian graphs shown in Figure |3] 
and Figure HI and the non-abelian graphs involving a three-gluon vertex shown in Figure 

We first discuss the abelian graphs. The two diagrams Dg and in the first row of 
Figure [3] are obviously a convolution product of two NLO functions and introduce no further 
computational complications. The sum of these two diagrams gives a symmetric color structure 



w- 



(8) _ _ 
ijk d 



{T?, If} T? Tl 



(35) 



and the integral is 



Ig(u, Qij, a ik ) = I [dk] [dl] 



Hi ■ rij rii ■ rik 5(u — n ■ (k + I)) 
rii ■ k rii ■ I rij ■ k n k ■ I 



2— 2e -2e^p -l-4c 

7T e 1E tu 



(36) 



with 



I 8 (a, a') 



g e 2e lE r 2(_ e ) r( ^_ 2e ) 

T(l -2e)r(l -4e) 
x 2 *i 



(1 -a)-' {I -a')~ e 

(-e, -e, 1 -e,a) 2 Fi(-e, -e, 1 



(37) 
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Figure 4: Example of a pair of mixed virtual-real one-particle cuts which adds up to a scaleless 
integral. 




(a) (b) 



Figure 5: Examples of non-abelian three- Wilson-line integrals required in the calculation of 
the NNLO soft matrix. 



Each of the four diagrams in the last two rows of Figure [31 on the other hand, are complicated 
functions of two distinct scalar products. However, the sums of the pairs (c) + (d) and (e) + (/) 
are proportional to symmetric color structure w^: 

(c)+(d) ^ (rrbrra , rrarrbx rrarrb f uu un n t -n 3 n t -n k 5(u — n ■ (k + I)) 



rii ■ I rii ■ (k + I) rij ■ k Uk ■ I 



D (e) + (/) ^ (T a rpb + T , ^ ^ f [dk] [dl] Vn 3 vn^ -V (H |)) 
8 v 4 * ; J J rii ■ k rii ■ (k + I) rij ■ k n k ■ I K ' 

Furthermore, after partial fractioning, the sum of the two integrals yields the factorized integral 
fl36l) . Therefore, these abelian diagrams do not introduce any new calculational complications. 
In Appendix [B] we explain how the non-abelian exponentiation theorem implies the simple 
factorized form of the integral multiplying the symmetric color structure (13"5|) . The color 
matrices for the qq channel are 
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while those for the gg channel are 
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We now discuss one-particle cuts of three- Wilson-line abelian diagrams. These are first of 
all due to contributions of the type shown in the first row of Figure EJ but where the connection 
to parton k is moved to the other side of the cut. These diagrams f actor ize in an obvious way 
into a product of an NLO-type integral with a scaleless virtual correction and thus evaluate to 
zero. The one-particle cuts corresponding to the diagrams in the last two rows of that figure 
are more complicated. However, one can show that the sum of pairs of the type shown in 
Figure H] reduces to the same scaleless integral mentioned above. Therefore, these one-particle 
cuts do not contribute to the soft function. 

In addition to the abelian diagrams, we must also evaluate non-abelian ones of the type 
shown in Figure However, the sum of diagrams in the figure vanishes due to its color 
structure. Written explicitly, the two diagrams give 

D' 8 {a) I? I* I? J£, 



D 8 {b) ->■ -if abc 



T t a T^T^I' 8 , 



where 



[dk][dl] 



M V P 



[(2k + l) v gp P - (k + 2l)^g pu + (I - k) p g ull ] 



71; 



k rij ■ I n k ■ (k + I) (k + iy 



(41) 
(42) 



(43) 
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Figure 6: Four- Wilson-line diagrams contributing to the NNLO soft function. 



The relative minus sign between the two diagrams originates from the fact that the gluon 
propagator carrying momentum k+l appears on opposite sides of the cut in the two diagram^. 
Therefore, the sum of the two diagrams cancels. This is due to the structure of the color 
factors and not the integral itself, and is also true of one-particle cuts not shown in the figure. 
Interestingly enough, it would also be true if the velocity vectors were time-like, even though 
such anti-symmetric three-particle correlations are present in the anomalous dimension for 
massive particles [4THI4T]. This is not a contradiction, since the soft function for massive 
heavy-quark production obeys an RG equation analogous to (j49j) below and the contribution 
of the three-particle correlations cancels between the sum of terms involving ~f s and 7]. 



3.3 Four- Wilson-line integrals 

We finally turn to the four- Wilson-line integrals. The two-particle cuts shown in Figure [6] are 
a convolution of NLO integrals and of the same form as I± and Is- The general integral reads 



Ig(uj, Qij, Qkl 



[dk] [dl] 



rii ■ rij rik ■ rii 5{oj — n ■ {k + /)) 
rii ■ k rij ■ k rik ■ I rii ■ I 



(44) 



and the color factor is 



d R 



(45) 



For two-to-two scattering considered here the relations between the different scalar products 
implies that we need only J 9 (o;, a^-, a^) = 7 4 (co>, a^), which we already took into account in 

- tuff in the qq channel, and 



(46) 



writing (1201) . The explicit results for the color matrices are 
for the gg they read 
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(9) 
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Ca 
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(n 2 -i 
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\ 



jV 2 -4 . 
' 47V 2 / 



4 This change of sign also occurs in the diagrams generated from the non-abelian two Wilson-line diagrams 
Dq and Dj 2 in Figure [2] by moving the two-gluon attachment to the other side of the cut. However, in those 
cases the order of the two color generators T" and Tf is also exchanged, which compensates the sign change 
of the propagators. 
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and w!?' = w$ for ij = 13, 14. 

Also in this case there are one-particle cuts, but these always involve a scaleless loop graph 
and vanish. 



4 Renormalization 

The renormalized soft function is obtained from the bare one by multiplying on both sides 
by a matrix-valued UV renormalization factor. The structure of this renormalization factor 
follows from the RG equation for the soft function, which is simplest to discuss at the level of 
the Laplace-transformed function. We define this as 



a (L, h/M 2 , fi) = jf duj exp {- — -^ j S(u, h/M 2 , fi) 
The integral transform is easily carried out from the bare function using 



(47) 



duexp[— bu]u' 



-1— lie 



T(-ne)b e 



(4c 



The RG equation for the Laplace-transformed function was derived in [12], using RG invariance 
of the cross section and the results for the hard, fragmentation, and parton luminosities. It 
reads 
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s ( In— ,ti/M 2 ,/i 
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A(a s )\n— + ~f s (ti/M 2 ,a s ) 



(49) 



where A = 2Cpj cvlsp in the qq channel and A = (Ca + CV)7cus P in the gg channel, and 

7 s (ti/M 2 , a.) = J h (h/M 2 , a s ) + [2 7 *(a.) + 2^«(a s )] 1 . (50) 

The non-trivial matrix structure of the soft anomalous dimension is related to that in the hard 
function and is expressed through the function 7 , which for the different channels is 



7^ = 47^)1 + ^7, 



cusp V^s J 
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(51) 



and 



7* = [27 9 («J + 2 1 "(a s )] 1 
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Explicit results for the expansion coefficients defined as 



(53) 



and analogously for the other functions can be found in the the Appendix of [T2] . 

Given the structure of the RG equations, we define the relationship between the bare and 
renormalized soft functions according to 



(54) 



The bare soft function does not depend on /x. This implies an RG equation for the renormal- 
ization factor which can be integrated to solve for its explicit expansion in e (see, e.g., [ST]). 
This yields 
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lnZ s = ^ + 
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(55) 



Defining expansion coefficients of the renormalization factor and Laplace-transformed func- 
tions in units of a s /ATc as in ([5]), the renormalized NNLO function in Laplace space is given 
by 



+ Z 



t(i)sW +s (1) 



bare 



bare 



Z (i) + z td) 5 (0) z (i) 



(i) 

'bare ' 



(56) 



Evaluating this equation, we find that the renormalized function on the left-hand side is indeed 
finite in the limit e — > 0. This shows that the renormalization factor ( 15 5|) following from the 
RG equation deduced from the factorization formula ([T|) is correct, or can otherwise be viewed 
as a cross-check on our calculations. The explicit expressions for the renormalized functions 
are rather lengthy and we relegate them to Appendix O The terms proportional to powers 
of L in the renormalized function are in agreement with the approximate NNLO formulas 
mentioned in [12], while the L-independent pieces are new. 
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5 Conclusions 



We evaluated the NNLO corrections to the soft function needed to describe the pair invariant 
mass distribution in top-quark pair production at hadron colliders in the small-mass limit. 
At the technical level, this required us to obtain real-emission type corrections related to a 
product of Wilson-loop operators depending on four light-like Wilson lines. This is the first 
NNLO calculation of a soft function which involves non-trivial matrix structure in color space. 
We showed that the IR structure of the bare function is consistent with known expressions 
for the two-loop anomalous dimension matrix derived in [T2]. The final results, given in 
the Appendix, turned out to be rather simple. This is because in the sum of all diagrams 
contributions from integrals involving three or more Wilson lines multiply a color structure 
whose coefficient is constrained by the non-abelian exponentiation theorem to be a product 
of NLO integrals. The non-abelian three-parton graphs are not constrained by non-abelian 
exponentiation, but these vanish after summing over all diagrams. 

Combined with known NNLO results for heavy-quark fragmentation functions and virtual 
corrections to two-to-two processes, our calculations will allow an evaluation of the invariant 
mass distribution in the small-mass limit at the level of a full virtual plus soft approximation. 
This will provide valuable information concerning the importance of higher-order corrections 
to NLO+NNLL predictions from [I]. Moreover, by an appropriate modification of the delta- 
function constraint in the definition 01]), we can immediately obtain the soft function needed 
to study single-particle distributions in the px and rapidity of the top quark to the same 
level of accuracy, opening up the possibility of studying higher-order corrections to the results 
of (121113]. Finally, we anticipate applications for threshold resummation in dijet production. 
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A Calculation of 1$ as an example 

Written out explicitly, the integral I5 reads (with 5 + (k 2 ) = 5(k 2 )6(k )) 



h(co,a tl ) = / d%d% 



(m ■ n 3 f 6+ (ll) 5 + (lp 6(00 - n ■ (h + Z 2 )) 
Ui -IxTii- (h + l 2 ) rij ■ lx rij ■ (7i + l 2 ) 
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d Uj ■ rij 5(lj -n -k) x _ t r(-e) ( 2 rii ■ k rij ■ k 



rii ■ k rij ■ k T(l — 2e) \ rii ■ n 

x (*) (-e,-e,l- e, 9 Wi 7*% ) • (57) 
v ' \ 2 n,i ■ k rij ■ k J 

To derive the equality in the third line we used the auxiliary integral given in [2^126]. The 
remaining d- dimensional integral over k can be solved using the steps outlined in |26j. We 
first define light-cone coordinates 

— l. i 1 U J 1 

ft — K + -f -f Hj_ , 

a/2 rii ■ rij a/2 rii ■ rij 

Tlj ■ k Tlj ■ k 2 _ 7,2 

y/rii ■ rij/2 y/rii • n.,/2 

and use these to parameterize the integral over k. This yields 

2 2 ~ 2e 7r 1 ~ 2e rY-f N )ri'i - e) r r°° r l 

I B (w,aij) = 1 ; 1 - \ dk+dk. / dk T k l T - 2e I dcos9 shr 1 " 2 ^ 

x (k + k_)- 1+e (k + k_ - k^)- 1 ' 2 " 2 F X ( -e,-e,l- e, 



( k + n _ + k_n 0+ \ 
x il u h n Tkr cos 6U . (5 c 

We then make the change of variables 

to to to 2 
k + = x y , k_ = (1 — x) y , dk + dk^ = y dxdy , 

kx = — \/ x(l — x) y u , dkx = V x (^- ~ x ) V d- u ? 

and express the light-cone components of uq in terms of to arrive at 

o2+2e 1-26 p( — ,\ p/i _ \ 

T(,,n \ - -S. rZ_A_ V - ,-1-46 /1 _ \-2e 
i 5 (U;, OyJ - — LU [l- Qij 



rl poo pi pi 

x / dx I dy dcosO sm~ l - 2e e / du[x{l - x)]' 1 - 26 y" 1 '^ u 1 ' 26 {I - u 2 ) 
Jo Jo J-i Jo 



2\-l-2e 



x 2 -Fi(— e, — e, 1 — e,u 2 ) 8\l — ^ + ^ a^xil — x)yu cos 6 J . (59) 

We next perform the y-integration using the delta-function constraint and integrate over cos 9, 
generating a hypergeometric function written in the first equality below. After that we perform 
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a change of variables to the argument of this hypergeometric function, integrate, and perform 
another change of variables to the argument of the original hypergeometric function. Explicitly, 



I B (uj,a 
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Jo 

The integral over the hypergeometric functions is sufficiently complicated that we evaluate it 
as an expansion in e. This is achieved using the standard identity 



oo 

i-r 1 - 6 =--*(i-*)+£ 



-ne) 



m=0 



ml 



ln m (l -t) 



(61) 



and evaluating the integrals over the plus distributions after expanding the other parts of the 
integral in e. This leads to parametric integrals over HPLs. Some of these are straightforward, 
and even the most difficult ones are not hard to handle by deriving and solving differential 
equations with respect to the parameter a^, yielding the result given in Section [31 



B Constraints from non-abelian exponentiation 

In this section we briefly explain the consistency of our results with the non-abelian exponen- 
tiation theorem j2ZH2B]- In general, this theorem states that the soft function can be written 
as the exponential of a simpler quantity. In particular, in the language of [28], the exponent 
receives contributions only from diagrams involving single connected webs. For us, the impor- 
tant point is that this theorem implies that after summing over all diagrams the coefficients of 
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certain color factors in the NNLO function are determined by exponentiating the NLO result. 
This exponentiation occurs in position or Laplace space rather than momentum space. We 
find it convenient to work with the Laplace-space function (I4"T|) . Our procedure is to take the 
exponential of the one-loop bare function and expand it to second order: 



1 



IR 



exp 
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d R 



+ 77 ^ T i ■ T j T k ■ Ti h(L, ciij) h(L, a H ) + 



i,j,k,l 



(62) 



The term on the second line produces two, three and four parton terms whose structure must 
be reproduced in our explicit results. First, we can write the two parton terms as 



T a . T a T 6 . T ^( L? ai2 ) + T « . ■ T b J({L, a 12 ) 



+2T? • T 3 a T 1 " • T*I 2 X {L, a 13 ) + 2T? ■ T?T? ■ T^l'f(L, a u 



b j2i 



(63) 



(3) 

We thus expect the coefficient of w\- to be restricted by the non-abelian exponentiation 
theorem, and a relation between J3, J4, I5 and the NLO integral. Indeed, we have already 
made explicit in (1221) that ^3 = 14/2-/5, and one can easily verify that I4 = If. After 
rewriting the color factor as in the last equality in ( 132]) and performing the sum over legs 
one can then show that the coefficient of w\j in ( l2~Tj) is proportional to the factorized integral 
I4, as required. The remaining two- Wilson- line integrals are proportional to the NLO matrix 
and are thus single connected webs which contribute to the second-order expansion of 
the exponent directly. 

Next we deal with three-parton terms. To understand their structure we consider the 
concrete example where two gluons attach to the parton with velocity rii and the other two 
to the partons with velocity n 2 and n^. The contribution of these terms to (1621) is 

2 ^ jiarpb . T b + T a . T a T b . j*j ai2 )j l(L; ai3 ) = 2w g/ 1 (L, a 12 )h(L, a 13 ), (64) 

where we have used that the color matrices for different partons commute. The above relation 
implies that h(L, ai2)Ii(L, 013) = Ig(L, ai 2 , ais), and one can confirm that this is indeed 
the case. Including all permutations in the expansion of the exponential, we reproduce the 
contributions proportional to the wfjl in (j2"Tl) . 

Finally, we turn to the four-parton terms, considering as an example the case where the 
two partons with velocity n\ and n 2 are connected to each other. These contribute to fl62l) as 

4 -Tf • T 2 a T 3 6 • T 4 b / 1 (a 12 )/ 1 (a 3 4) = Aw^h(a 12 ) . (65) 



d 



R 



This accounts explicitly for the four-parton term wf^ I 4(^12) (|2"T|) . The remaining four-parton 
terms appear as wff J 4 (ay). 
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C Renormalized Soft Functions 



We list here the results for the renormalized soft function in the qq and gg production channels. 
For the sake of brevity we set iV = 3 and take into account that the soft function is symmetric. 



In the following, we indicate the element ij of the matrix s^f 1 (k G {qq, gg}) as sj^- 



s(») 



C.l Quark Annihilation Channel 

The elements of the NLO soft matrix in Laplace space are 
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The elements of the NNLO soft matrix in Laplace space are 
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1 L H iV; 

3 3 9 



rr / \ 160 rr . , 128 rr , , 160 rr . , 
i^ifo) - — H 2 ,2(x t ) + -y-H 3fi {x t ) + —H 3>1 (x t ) 



+ 



176 448 r 32 Ar 

+ L + —Ni 

3 3 9 



H 0> o,o(x t ) + 



176 352 r 32 Ar 

+ L + —N 

3 3 9 



Hi, 0>0 {x t ) 



+ 



176 160 r 32 

+ L + — iVj 

3 9 9 



Hi,i,o(xt)+ 



176 128 r 32 Ar 

1 L H Ni 

3 3 9 



448 TT . . 160 rr . . 416 rr . . 736 TT . . 160 rr . . 
- —^i,i,2(^) - -g-tfiAofct) - — ^1,2,1 (a: t ) + — #2,0,0(2*) - —H 2 ,i,o(x t ) 



64 rr . . 448 rr . , 256 rr , , 160 rr , , 
+ y#2,i,i(z*) + —H ,oM x t) + —H h0fi ,o{x t ) + —H hh0fi {x t ) 



23 



256 128 
- — #i,i,i,o(zt) - — #1,1,1,1 ( x t) , 

(2) 38848 12928 2144 2 352 3 256 4 5248 !792 

s W 22 --243 ^r L + ^r L ~Yf L + Yf L ~T2a Nl + ^m LNl 

320 .2 ^ , 64 268 2 308 2 416 40 , 56 , 



994 
405 



vr 4 + 



11312 3 75 2 308 r2 896 , 1568 Ar 560 r „ 56 r2AT 

1 L L H L H N t LNi H L N, 

81 27 9 27 243 81 27 



308 2 856 r 2 56 Ar 2 392 A 
27* + ^ + 8T^ + - Cs 



#o(&t) + 



3232 1072 



81 



27 9 27 



448 Ar 160 rAr 16 r2Ar 308 2 80 r 2 56 Ar 2 

JVi H LiV/ L Afj 7T 2 Lvr 2 H N^ 2 

243 81 27 27 81 81 



#i(z t ) + 



1072 



176, 64 , 160 .. 32 



-L + -L 



80 



81 



Ni + —LNi + — 7T 



27 



81 



176 64 

1 L 

9 9 



32 
27 



Ni 



H 3 (x t ) 



32 



3752 616 r 2720 

L H , 

27 9 27 



560 



112 



856 



N t H LNi H 7T" 

81 27 81 



H ,o(x t ) 



3752 616 , 704 



560 



112 



352 



27 



-L + L N t + LNi + 

9 27 81 27 81 



7T 



Hi, (x t ) 



1072 176 



27 9 



640 
~27 



160 
~8T 



L JV, H L/V; 7T 2 

27 



424 
~8T 



Hi,i(x t ) + 



176 



128 r 32 Ar 

L Ni 

9 27 



64 

+ y#i, 3 (^) + 



616 608 r 112 Ar 

1 L H Ni 

9 9 27 



#2,0 (x t ) + 



176 64 r 32 

L jVj 

9 3 27 



#2,l(^) 



160 TT . . 736 rr . . 128 rr . . 
+ —#2,2(2^) + — iJ 3 ,o(^) + —H 3<1 (x t ) + 



616 608 112 Ar 

1 L H Ni 

9 3 27 



#o,o,o (xt) 



+ 



616 _ or 112 Ar 
+ 128L + —N t 



27 



616 64 r 112 Ar 



Hi,i, (x t ) 



+ 



176 128 r 32 Ar 

L Ni 

9 3 27 



rr / \ 256 rr , , 400 rr , , 320 TT , , 
#i,i,i(^t) + — #1,1,2(3*) + -g-#i, 2 ,o(zt) + -g-^i, 2 ,i(^) 
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1280 TT , v 64 TT . . 64 rr , , 608 TT . N 

+ -g--H2,0,0(Zt) + y# 2 ,i, (Zt) + y #2,1,1 0&t) + — H OfiM X t) + 128^1,0,0,01^) 

608 rr , . , . 128 TT , . 352 448 r , 64, r> 

+ -g-#i,i,o,o(z t ) " 16#i,i,i,o(z t ) + -g-ffi,!,!,!^*) - y^3 + — ^3 + . (67) 

C.2 Gluon Fusion Channel 

The elements of the NLO soft matrix in Laplace space are 

S£> u = 26L 2 + % 2 , 

& = 2 * 2 + l2LH (x t ) + ULH^Xt) - 12H 2 (x t ) + 12H 0fi (x t ) + 12H lfi (x t ) - l2H hl (x t ) , 

s gg 13 - u > 

4s22 = 13L ' + + 9^o(^) - QLH^Xt) + 9if 2 (x t ) + 9# ,ofa t ) + 9H lfi (x t ) + 9H 1A (x t ) , 

Sgg 23 = l* 2 + 5LH (x t ) + hLH^Xt) - 5H 2 (x t ) + 5H 0>0 (x t ) + 5H lfi (x t ) - 5H 1;1 (x t ) , 
55 65L 2 

s 9 g33 = ^ 2 + -jp + 5LH (x t ) - 5L# 1 (x t ) + 5# 2 (x t ) + 5H 0fi (x t ) + 5H 1>0 (x t ) + 5H 1A (x t ) . 

(68) 

The elements of the NNLO soft matrix in Laplace space are 



(2) 31564 10504 1742 2 286 r3 338 r4 4264 1456 AT 260 
* M n = -27 9~ L + ^~ L + ~ L " "8T^ + "27^ " ~ J 



52 o 871 9 104 , 9 65 , 461 4 9 , N , 

+ yL 3 ^ + — TT 2 + y-vr 2 L 2 - -it 2 N x - — vr 4 + 96L 2 H 2 (x t ) + 96L 2 H 0fi (x t ) 

+ 9QL 2 H h0 (x t ) + 9QL 2 H hl (x t ) + 167r 2 L# fat) + 16rr 2 Lifi(a; t ) - 96LH 3 (x t ) 

- 96LH 1:2 (x t ) + 96LH 2 , (x t ) - 288LH 2A (x t ) + 288L# ,o,ofat) + 288LH 1:0 , (x t ) 

+ 96LH hlfi {x t ) - 288LH w {x t ) - 16n 2 H 2 {x t ) - 288H 4 {x t ) + 167r 2 # ,ofat) 

+ 167r 2 // 1>0 (^) - 167r 2 // 1;1 (x t ) - 288H h3 (x t ) - 96H 2:2 (x t ) - 9QH 3fi (x t ) 

+ 96# 3 ,ifat) - 96^i,i,2 (x t ) - 96// 1>2)0 (x t ) + 96^i )2 ,i(x t ) + 96# 2 , ,ofa t ) - 288if 2 ,i,o(a;*) 

+ 2S8H 2tl>1 (x t ) + 288# ,o,o,ofat) + 288# 1 , OA ofa t ) + 96^i,i )0l o(a; t ) - 288/J 1 , 1 , 1 , (^) 

S2 286 
+ 288# 1>1 , 1 , 1 (z t ) + 1092LC(3) + -N,C(3) - — C(3) , 
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r(2) 
'9912 



134 „ 52 



4 
3' 



-TT 2 + — 7T 2 L 2 + -TT 2 LNi - 227T 2 L - — 7T 2 N t - —TT 4 



20 



6£ 
45' 



104L 3 + AL 2 Ni - QQL 2 



fLN t + ^rr 2 L + 268L + ^ + H?JV, + 252C(3) - 22tt 2 - ^ 



#o(**) 



+ 



40 16 112 4 

104L 3 + AL 2 Ni - Q6L 2 - —LN t + —n 2 L + 268L + — N t + -vr 2 ^ + 324C(3) 

3 3 9 3 



- 22tt 2 - 



308 



40 



16 



-10AL Z - 8LNi + 132L + —N t n 2 - 268 



H 2 (x t ) 



+ 



40 



52 



176L Z + 8LNi - 132L - — JV, + — tt z + 268 



40 16 9 

y iV, + yTT 2 + 268 



Hifl(xt) 



104L 2 + 8LJV, - 132L 



40 „ 52 



-176L" - 8LiV, + 132L + — iV ; - —n z 



268 



H hl (x t ) + [72L + 8^ - 132] H lj2 (x t ) + [72L + 8iV, - 132] H 2fl (x t ) + 



1AAL 



+ 8iY, - 132 



H 2 ,i(x t ) + 



216L + 8N, - 132 



+ 



216L + 8Ni - 132 



H 1X i{x t ) + 



H 0fi ,o(x t ) + 



H 3 (x t ) + 



1AAL + 8N1 - 132 



Hifi,o(x t ) 



8N[ - 132 



Hi,i, (x t ) 



8N t - 132 

+ 72H 4 (x t ) + 72H 3fi (x t ) - 72H^ 2 {x t ) - 144if 1)2)1 (a; t ) + lAAH 2fifi (x t ) 
- 72H 2:hl (x t ) + 216# ,o,o,o(zt) + 72H lflfifl {x t ) - 72H Wfi (x t ) - 21QH 1AAA (x t ) , 
(2) 5tt 4 20 



H (x t ) + H^xt) 



r2 20 2 
40L 2 - — 7T 2 

O 



H 2 (x t ) + H^xt) 



+ 



40L 2 + ^tt 2 



#o,o(^) + H h0 {x t ) - 40L# 3 (* t ) - 40L#i, 2 (z t ) + A0LH 2fi (x t ) 



- 120Lff 2)1 (x t ) + 120L# ,o,o(z t ) + 120L// 1AO (x t ) + 40Lif 1)li0 (a; t ) - 120L#i,i,i(:r t ) 

- 120# 4 (x t ) - 120if 1)3 (a; t ) - 40# 2 , 2 (z t ) - 40# 3 , (a; t ) + A0H 3A (x t ) - A0H 1A:2 (x t ) 

- A0H lt2fi (x t ) + A0H 1;2A (x t ) + 40// 2 ,o,o(^) - l20H 2>lfi (x t ) + 120H 2 , 1;1 (x t ) 

+ 120H OA o t0 (x t ) + 120^,0,0,0^) + 40^1,1,0,0(0:*) - 120^1,1,1,0^) + l20H WA (x t ) , 



26 



(2 ) 15782 5252 T 871 2 143 rS 169 r4 2132 Ar 728 Ar T 130 Arr2 
' ~ L + —L 2 - —L 3 + —L A - -^-Ni + NiL N t L 



>gg 22 ~ 2 7 

26 



9 

2077 



81 



27 



-7T 



33 o T 169 o ^ o 

— iTL + —7rL 2 



155 
~54 



■iV/7T 2 + jVjTTL 



673 
135 



9 



7T 4 + 



78L d 



99 33 28 

+ 3iV,L 2 L 2 - 10JV,L + 26tt 2 L + 201L + JV,7r 2 vr 2 - 202 + — N t 

2 2 3 



+ 189C(3) 



H (x t ) + 



-78L 3 - 37V,L 2 + ^L 2 + WNiL - 201L - — JV, + jV,7r 2 



33 

- — vr 2 + 202 - 135C(3) 



H 1 (x t ) + 



102L 2 + 6JV,L - 99L - 10Ni + 201 



+ 



210L 2 + 6N t L - 99L - 10iV, + 26tt 2 + 201 



#o,o (&t) + 



102L 2 + 6N t L - 99L 



- lOiV, + 17tt 2 + 201 
+ [-24L - 6JVj + 99 



#i,o(z t ) + 



210L 2 + 6N t L - 99L - 10N t - 9vr 2 + 201 



H 1A (x t ) 



H 3 (x t ) + [-132L - 6iV, + 99] #i, 2 (^) + [-288L - 6N t 
+ 99] # 2 ,i(x t ) + [l32L + 6iV, - 99] # 2 , (x t ) + [396L + 6iV ; - 99] # ,o,o(^) 
+ [288L + QNt - 99] H lfi , (x t ) + [24L + 6iV z - 99] H hlfi (x t ) + [-396L - QN t 



+ 99 



#i,i,i(z t ) - 126H 4 (x t ) - 72H h3 (x t ) + 30H 2 , 2 (x t ) + 84H 3fi (x t ) + 132H 3>1 (x t ) 



+ MH 1X2 {x t ) + 30#i, 2 , (z t ) + 2A0H lt2tl (x t ) + 240H 2fii0 (x t ) - 72H 2tlfi {x t ) + 2S8H 2X1 (x t ) 
+ 396H mo (x t ) + 2S8H lfit0t0 (x t ) + 132if lil)0l0 (a; t ) - mH lxlfi (x t ) + 396H ltlxl (x t ) 

-^C(3) + fiV,C(3) + 546C(3)L, 



4?23 = ^ 2 ~ + ^ 2 L 2 - ^N t rr 2 + "-N^L - ^tt 4 



335 



55 
~6~ 



65 



25 
27 ^ 



5 
9 J 



23 
108' 



130 



55 

y 



Arr 110 2r 335 r 140 Ar 5 Ar 2 55 2 1010 inr . /0 , 
A^L + — vr 2 L + — L + — iV, + -JV,tt 2 - -vr 2 - — + 105C(3) 



H (x t ) 



+ 



130, 5 Arr2 55 r2 50 Arr 20 2 r 335 r 140 5 Ar 2 55 2 1010 

L + -NiL 2 L NiL H n 2 L H L H iVj + -Niir 2 n 2 

3 3 2 9 9 3 27 9 6 9 



27 



+ 135C(3) 



H 1 (x t ) + 



130 r2 10 Arr rrr 50 Ar 20 2 335 

L NiL + 55L H N t tt 2 

3 3 9 9 3 



H 2 (x t ) 



310 



10 Arr rrr 110 2 335 50 Ar 
+ yAT / L-55L + — vr 2 + — - -JV, 



H ,o(x t ) + 



130 o 10 50 

L - 55L + — NiL Ni 

3 3 9 



65 2 335 



Hi, (x t ) 



310 r2 10 Arr rrr 50 Ar 65 2 335 

L NiL + 55L H N, tt 2 

3 3 ' 9 ' 9 3 



H 1A (x t ) 



+ 



60L + —N t - 55 
3 



H h2 (x t ) + H 2fi (x t ) 



120L + —Ni - 55 
3 



#2,1 (x t ) + 



180L 



+ y N t - 55 



H ,o,o(x t ) + #i,i,i(z t ) 



+ 



120L + — Ni - 55 



Hi >0 ,o(xt) + 



10 



55 



H 3 (x t ) + 30H 4 (x t ) - 30H 2 , 2 (x t ) + 60# 3 , (x t ) - 60H 3>1 (x t ) + ^-NiH hlfi (x t ) 



- 55H 1X0 (x t ) - mH lX2 {x t ) + 30ifi )2)0 (a; t ) - 120H h2<1 (x t ) + 120H 2fifi (x t ) 

- l20H 2>ltl (x t ) + 180# ,o,o,o(z*) + 120^,0,0,0^) + 60ffi i i, Ol o(a; t ) - 30H 1AAfi (x t ) 

- lB0H 1>ltltl (x t ) . 



(2) 78910 26260 T 4355 r2 715 r , 845 r4 10660 Ar 3640 Ar T 650 Arr2 
( ' ~ L + —r^L 2 - —L 3 + —L 4 - ^r—Ni + N t L N t L 2 



99 33 243 

130 



27 



27 



27 



729 



243 



81 



10385 2 55 

7T 

324 6 



845 
~8T 



+ N t L d + -tt 2 7T 2 L + 71 L - —Nili 1 + -Ni7l Z L - — 71* 



775 
486 J 



5 
9 J 



763 



243 



130 



5 Arr2 55 r2 335 T 50 Arr in 2r 5 Ar 2 55 2 1010 140 Ar 

+ -N t L L + L N t L + 10tt 2 L + -iV/Tr 2 tt 2 + — N t 

3 2 3 9 9 6 9 27 



+ 105C(3) 



H (x t ) + 



130 r3 5 Arr2 55 r2 50 Arr 40 2r 335 r 5 Ar 2 

L 3 - -N t L H L H iV z L n 2 L L + -iV/vr 2 

3 3 2 9 9 3 9 



55 

~6~ - 



1010 140 



9 -^-^(3) 



H^xt) + 



30L 2 + —N t L - 55L - -iV, + -tt 2 
3 9 9 



+ 



335 



H 2 (x t ) + 



90L 2 + ^NL - 55L - ^N t + 10vr 2 + ^ 
o y o 



#o,o (x t ) + 



30L 2 



28 



10 Arr rrr 50 Ar r 2 335 

H NiL - 55L Ni + 5ir H 

3 9 3 



+ 



335 



H 1A (xt) + 



H h o(x t ) + 



H 3 (x t ) 



90L 2 + -N^ - 55L - — JV, - -vr 2 
3 9 9 



140 r 10 Ar rr 



- #1,2(2*) 



- #1,1,1 (ac t ) 



10 



+ 



-80L - — JV, + 55 



30L + — N t - 55 
3 



H 2> i(x t ) 



#1,0,0(2*) + 



UOL + —Ni - 55 



#2,0(2*) 



#0,0,0(2*) 



-^L-55 



#1,1,0(2*) + 10H 4 {x t ) 



+ 40#i, 3 (z t ) + ^# 2 , 2 (a*) + ^#3,0(2*) + ^#3,1(2*) + yiVjtfi.i.ofo) 

220 130 320 320 

+ — #i,i, 2 (x t ) + — #i, 2)0 (a*) + — #i, 2 ,i(x t ) + —#2,0,0(2*) + 40# 2 ,i,o(x t ) 

140 

+ 80# 2 ,i,i(zt) + 140# , ,o,o(2*) + 80#i, ,o,o(^) + —#1,1,0,0(2*) + 10#i,i,i,o(z t ) 

f v 130^C(3) 715 > , . 910L , . 
+ WOFm,!^) + — - — C(3) + — C(3) • 



(69) 
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